Spin currents in thermodynamic equilibrium? 
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The standard definition of a spin current, applied to the conductors lacking inversion symmetry, 
results in nonzero spin currents. I demonstrate that the spin currents do not vanish even in ther- 
modynamic equilibrium, in the absence of external fields. These currents are dissipationless and 
are not associated with real spin transport. The result should be taken as a warning indicating 
problems inherent in the theory of transport spin currents driven by external fields. 
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Generating spin currents (SCs) is one of the central 
goals of spintronicsfi and various mechanisms for elec- 
trical and optical injection of SCs have been proposed. 
Recently, the interest in SCs has been strengthened by 
independent reports of dissipationless SCs in two differ- 
ent systems: holes in the valence band of a diamond type 
crystal described by a Luttinger Hamiltonian,- and elec- 
trons in a 2D system with a structure inversion asym- 
metry (SIA) described by a Rashba Hamiltoniani^ These 
systems differ in symmetry because the Luttinger Hamil- 
tonian possesses inversion symmetry while the Rashba 
Hamiltonian lacks it. In both cases, SCs are driven by an 
external electric field E . These surprising results have 
drawn a lot of attention, causing active interest and im- 
mediate response. They were followed by several papers 
of different researchergii^ and also by the more recent 
papers coming from the same groupsiSiS The mathe- 
matical formalism in some of these papers, particularly 
in Refs. 2a 5 and 2o is rather involved. Meanwhile, the 
somewhat miraculous nature of the dissipationless SCs 
calls for a better understanding of their mechanism, in- 
cluding the properties of the background that supports 
the SCs linear in E . From the standpoint of spintronic 
applications, it is important to understand whether these 
are transport cwcTents, i.e., whether they can be employed 
for transporting spins, accumulating them at specific lo- 
cations, and for injecting spins. 

In this article I am trying to contribute to this basic 
physical understanding. I consider non-centrosymmetric 
2D systems in thermodynamic equilibrium and show that 
using the standard definition of a SC results in non- 
vanishing SC expectation value. This result is not en- 
tirely surprising from a general symmetry viewpoint. In- 
deed, the reversal of the momenta (or the velocities) only, 
without reversing the angular momenta (spins), cannot 
be reconciled with the time-inversion symmetry for non- 
centrosymmetric systemstiS Of course, such background 
currents, which are present in the ground state in the ab- 
sence of in-plane external fields, are non-transport cur- 
rents. 

The standard Hamiltonian with a Rashba term is 

Hr = h'^P/2m + aR{cr x fc ) • z, (1) 
with cr the Pauli matrices, k ~ {kx,ky) the 2D mo- 



mentum, and z a unit vector perpendicular to the con- 
finement plane. The eigenvalues of the Hamiltonian Hn 
are E\{k) = h^k"^ /2m + X\aR\k , where A = ±1 corre- 
spond to the upper and lower branches of the spectrum, 
respectively. The eigenfunctions are 



"^^^ ^ " V2 ( -i>^(^Rik. +iky)/\aR\k) ) 
The operator of the velocity is 

V =h~^dHR/dk = hk /m + aji{z X cr), 



(2) 



(3) 



and the mean values of the Pauli matrices in the cigcn- 
states •ip\{k) are 

(,T),fc ^{^^[k)\cT\Mk))^^{k xz). (4) 

These equations allow one to evaluate the SC compo- 
nents defined as products of the electron velocity and 
spin components. 

For calculating SCs in a given state (A, A; ) , it is con- 
venient to evaluate the Hermitian dot and cross products 
of cr and v by using Eqs. © and 10}: 

]^[{a ■ v) + {v CT)]^^{a ■ k), {{a ■ k))^j^ ^Q, 
and 



\{[{crxv)-{v X ^)]),fc = f (2 + \^\ 
2 Ti \ ■m\aR\J 



(5) 



The dot-product vanishes because the spin is polarized 
perpendicularly to k in the eigenstates V'a ( ^ ) ■ The 
Kramers conjugate states, ip\{k) and %p\{—k), be- 
long to the same branch of the spectrum. As a re- 
sult, SCs of Eq. (O are even with respect to the change 
k — k within each branch, while their values for 
different branches, A = ±1 , are not mutually related. 
Therefore, there is no compensation of the SCs from dif- 
ferent branches, and net macroscopic SCs can arise. 

Electron spin is not conserved in the presence of a spin- 
orbit interaction because of precession in a momentum- 
dependent effective magnetic field. Therefore, the usual 
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procedure for deriving currents from the continuity equa- 
tions is not applicable to SCs (at least in its simplest 
form), and I use in what follows the standard and phys- 
ically appealing definition of the SC tensor J7ij 




Here i,j — x,y , with i indicating the spin component 
and j the transport direction. For T = , the integra- 
tion should be performed over k < K± , where K± are 
the Fermi momenta for both spectral branches. For the 
Hamiltonian Hn the tensor J7ij is antisymmetric 

*y^xx ^yy t^^xy i^ya: — t^^R- ('^) 

This tensor is invariant under the operations of the sym- 
metry group C ooD of the Hamiltonian Hn 

When the electrochemical potential /i is positive, /x > 
, both spectral branches are populated and 

Jr{p) =m'^ay3nh^. (8) 

Therefore, the SC is odd in the coupling constant aji , it 
is of the third order in an , and does not depend on ^ . 
For comparison, the spin-orbit energy found by averaging 
the second term of Hr is even in an { fJ,> 0): 

EM - ~iA/3TT){maR/h^fima%/h^ + 3^/2). (9) 

For small electron concentrations, when /x < and 
electrons populate only the lower branch of the spectrum. 

Equations (jHJ and match smoothly at ^ — 0; a dis- 
continuity exists only in the second derivative. Near the 
bottom of the spectrum, at fi = i^min = —maj^/2h^ , 
•JxyifJ-) shows a square-root singularity. The nonana- 
lytical behavior at these points emerges because of the 
spectrum singularities at E — and E = i?min and is 
also known for different phenomena. 

The SCs of Eqs. ||SJ| and H1U|) were found under the con- 
ditions of thermodynamic equilibrium. Therefore, they 
do not describe any real transport of electron spins and 
cannot result in spin injection or accumulation. Calcu- 
lating transport currents would require a modification of 



Eq. ©. Nevertheless, Eq. Q provides some insight into 
the effect of spin-orbit interaction on an equilibrium elec- 
tron system. A real (i.e., invariant under time inversion) 
antisymmetric pseudo-tensor J7ij is equivalent to a real 
vector P 1 1 z . It has the symmetry of the normal elec- 
tric field E ± producing the SIA and can be related to 
the spin-orbit contribution to the response of 2D elec- 
trons to this field. Interestingly, electric fields generated 
by SCs in magnetic insulators were discussed by Meier 
and Lossii^ (Despite a similarity, both the origin and the 
scale of the effect are very different from the above.) 

Calculating SCs driven by an in-plane field E is be- 
yond the scope of the present work. Nevertheless, we 
point out at a resemblance between the background SCs 
Jxy — —Jyx of Eq. 13) produced by the field E _l and 
the dissipationless SCs of Refs.|3 and|^ Applying a driv- 
ing field £^ to a diamond type crystal violates its inver- 
sion symmetry and, therefore, can produce background 
currents. To obtain transport SCs, these background cur- 
rents should be eliminated. 

The symmetry of the tensor Jij depends on the spe- 
cific choice of the spin-orbit interaction. When it orig- 
inates from the bulk inversion asymmetry (BIA), the 
Hamiltonian includes a Dresselhaus spin-orbit term. In 
the principal cubic axes it reads 

Hd = h^k'^ /2m + aDi<Jxkx — cTyky). (11) 

As distinct from Hr , Hu possesses only C 2v sym- 
metry. As a result, the ordering of electron spins with 
respect to k is quite different, and 

t!Txx — ^yy — i^Td 1 ^xy — ^yx — 0, (f^) 

where Jd(^) can be found from Jr{h) of Eqs. © and 
((Tn|l by the substitution ur ^ ■ Therefore, Eq. ® 
results in equilibrium SCs also for BIA systems. 

In conclusion, the standard procedure for calculat- 
ing spin-currents, when applied to non-centrosymmetric 
crystals, results in non- vanishing currents even under the 
conditions of thermodynamic equilibrium. For calculat- 
ing transport spin-currents, a procedure for eliminating 
the background currents should be devised. 
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